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ABSTRACT 

We report the Spitzer Infrared Spectrograph (IRS) observations of three 
evolved stars: IRAS 16456-3542, 18354-0638, and 23239+5754. The 9.9-37.2 /xm 
Spitzer/IRS high resolution spectra of these three sources exhibit rich sets of 
enstatite-dominated crystalline silicate emission features. IRAS 16456-3542 is ex- 
tremely rich in crystalline silicates, with > 90% of its silicate mass in crystalline 
form, the highest to date ever reported for crystalline silicate sources. 



Subject headings: stars: AGB and post-AGB 
infrared: stars 



stars: circumstellar matter 



Introduction 



The field of circumstellar silicate research started in the 1960s when iKamijol (jl963[ ) 
first theoretically investigated the condensation of dust grains in the circumstellar envelopes 
around M type long period variables and argued that silica (SiO^) should be the most 
abundant condensed species. iGilmanl ( 119691 ) showed several years later that grains around 
oxygen-rich cool giants are mainly silicates such as A^SiOs and Mg2SiO"4. The presence 
of silicate dust in evolved stars was first revealed throu gh the detection of the 9.7 jum Si-0 
stretching emission band in M stars (jWoolf fc Neylll969l ). This was later fur ther supported by 
the d etection of the 18 fim O-Si-0 bending emission band, also in M stars ( jTreffers fc Cohen 
1974J ). The 9.7 and 18 /im silicate bands have also been seen in absorption in evolved stars 
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with an extended, heavily obscured circumstellar dust shell (see lJones fc Merrilllll976l ; lBedijn 
19871 and references therein). 



All the earlier observations showed that the 9.7 and 18 /im circumstellar silicate bands 
were broad, smooth and featureless, suggesting that they were predominantly amorphous 
(IMolster fc Kemperi [io05h I ndeed, they can be we ll fitted with laboratory data for amor- 
phous olivine silicates (e.g. see lDorschner et~aL 1995 ). This picture has dramatically changed 
in the 1990s, thanks to the Infrared Space Observatory (ISO) and the Spitzer Spac e Tele - 
scope. Using the Short Wavelength Spectrometer (SWS) on board ISO, IWaters et al.l (119961 ) 
obtained the 12-45 /im emission spectra of six oxygen-rich evolved stars. They for the first 
time reported the detection of crystalline silicates in evolved stars, as revealed by the narrow 
emission features characteristic of crystalline silicate minerals at wavelengths between 30 and 
45 /im. 

Crystalline silicates have now been observed in all evolutionary stages of evolved stars: 
red giants and supergiants, asymptotic giant branch (AGB) sta rs, post-AGB stars , plan 



Svlvester et al. 


1999: 


Molster et al. 


2002a 


b 



been found in evolved stars by Spitzer in the L arge Magellanic Cloud (e.g. iKastner et al. 



2006; Sl oan et al.ll2006l . 120081 : IZiilstra et ali2006f ) and in R 71, a LBV star by ISO flVoors et al 



19991 ). We note tha t like amorphous silicates, crystalline silicates are usually found in oxygen- 
rich environments ( Molster fc Kemper 20051 ) PI Howev er, the presence of crys t alline silicates 



in ca rbon-rich evolved stars has also been reported (IWaters et al.l Il998al ibl; iMolster et al. 



20011 ). In these cases th e crystalline silicate grains were likely pr oduced in previous oxygen- 
rich mass-loss episodes (IMolster fc Kemperi 120051 : iHenningi 120101 ). 



The identification of circumstellar crystalline silicates has been made based on the char- 
acteristic mid-IR bands of crystalline silicates. The ISO SWS and LWS (Long Wavelength 



1 The atoms in an amorphous silicate grain are arranged in a random manner. In contrast, in crystalline 
silicates, the highly ordered counterparts of amorphous silicates, the atoms are arranged in an ordered lattice 
structure. 

2 LBV stars are very bright, massive, short-lived, hypergiant variable stars, with a luminosity exceeding 
10 6 Lq, a mass up to ~ 150 M@, and a lifetime of only a few million years. LBVs are a stage in the evolution 
of very massive stars. They may evolve to Wolf-Rayet stars before exploding into supernovae. 

3 In an oxygen-rich environment, the C/O ratio is smaller than one. All car bon atoms will be bound in CO 
and there will be some extra oxygen atoms left to form silicates and oxides (jGaill 120101 ) . In contrast, if the 
circumstellar envelope around an evolved star is carbon-rich (i.e. C/O > 1), all oxygen atoms will be bound 
in CO and therefor e no silica tes or oxides will be able to form; instead, one expects to form carbonaceous 
dust and SiC dust jGaillboiol) . 



3 



Spectrometer) spectra of evolved stars in the 2.4 to 195 /mi wavelength region allow one to 
identify at least 49 narrow bands which are attributed to crystallin e silicates, with disti nct 
emission bands at approximately 10, 18, 23, 28, 33, 40 and 60 /im jMolster et alJl2002al fbh . 
The widths and peak wavelengths of these bands and their relative strengths contain useful 
information about the crystalline silicate composition (e.g. the Fe/Mg ratio) and stoichiome 



try (e.g. olivine Mg 2x Feo -2 X Si04 or pyroxene Mg x Fei_ x Si03 where < x < 1; see lKoike et al. 
19931 ; IJaeger et al.lll998l ). The mid-IR spectra of evolved stars obtained by ISO and Spitzer 



so far indicate that circumstellar crystalline silicates all appear to be ext remely Mg-rich and 
Fe-poor or even Fe-f ree (i.e. forsterite Mg 2 Si04 or enstatite MgSiOs; 
2005l : lHennin3l2010h . 



sec 



Molster fc Kemper 



However, there are several outstanding problems about crystalline silicates remain to 
be addressed: First, the exact composition and stoichiometry of crystalline silicate dust in 
evolved stars are not precisely known. Second, the mechanism for making circumstellar 
crystalline silicates remains debated. Third, it is unclear whether (and how) the fraction 
of silicate dust in crystalline form depends on the mass loss rate, the temperature of the 



circumstellar dust shell, and the evolut i onary stage of an evolved star (IKemper et al.l 12001 



Molster fc Kemperll2005l ; IWoodenll2008l ; lHenning||2010l ). It has also been suggested that high 
crystallinity (> 10%) is typically associated with the presence of a disk, rathe r than being 
correlated with the mass- loss rate ( Kemper et al.ll200ll ; lMolster fc Kemperll2005l ). Further, it 
is widely believed that interstellar dust originates from Stardust produc ed by evolved stars 
It is puzzling that interstellar sil icates are predominantly amorphous (IKemper et ah||2004] ; 
Li &: Drainei l200ll ; iLi et al.l 120071 ) while an appreciable fraction of circumstellar silicates - 
the primary source of interstellar dust - are crystalline. 

To improve our understanding of these questions, we plan to systematically explore the 
spectroscopic properties of astronomical crystalline silicates and their formation, evolution, 
and dependence on the properties of the parent stars as well as their relevance to inter- 
stellar silicates. In this work we present the Spitzer/IRS spectra of three evolved objects 
(IRAS 16456-3542, 18354-0638, and 23239+5754) all of which display rich sets of crystalline 
silicate emission features. This is a part of our ongoing efforts for a systematic investigation 
of crystalline silicates in evolved stars. 



4 But also see brainel (2009) and IZhukovska et al.l (|2008l ) who argued that the dust in the interstellar 
medium (ISM) is not really "stardust" but condensed in the ISM as the inter stellar dust destruction rates 
are higher than the injection rates of Stardust into the ISM ([Jones et al.lll994l ). 
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2. Observation and Data Reduction 

The data were obtained through the Spitzer Cycle 3 GO proposal (Spitzer Program 
ID 30403; PI: Biwei Jiang) fl Sixteen evolved objects (post-AGB stars and PNe) have been 
observed through this program. They were selected from the Infrared Astronomical Satellite 
(IRAS) Point Source Catalog (PSC), with their IRAS color indices satisfying ^(-^25/^12) > 
—0.2 and lg(F 6 o/-F 2 5) < 0. This is to exclude young stel lar objects (YSOs) or galaxies 



and limit us to evolved stars (Ivan der Veen fc Habing]|l988l ). In this work we focus on the 



three sources (IRAS 16456-3542, 18354-0638, and 23239+5754) whose mid-IR spectra are 
characterized by rich sets of crystalline silicate emission features. The other (13) objects will 
be discussed elsewhere. In Table [T] we list their IR brightness (and thus colors) and their 
SIMBAD classification. 

These sources were observed using the Infrared Spectrograph (IRS) on board Spitzer 
with the Short-High (SH) and Long-High (LH) modules during Cycle 3 from June 2006 to 
May 2007. The data were taken in 'staring mode' at two nod positions along each IRS slit. 
Thanks to the anticipated high brightness of the sources, no specific off-set observations for 
measuring the background were needed. The journal of observation is shown in TableEJ 

We start the data processing from the Spitzer Science Center's pipeline data set (version 
S16.1.0). The irsclean program is used to combine unflatfielded image (droopres) and the 
rogue pixels image from the same campaign for a compound mask, afterwards this combined 
mask is applied to the flatfielded image (bed) to remove rogue and flagged pixels as much as 
possible. The spectra are extracted from the cleaned image in the SPICE software package 
and further analysis is done in the SMART (version 6.2.6) package. 

Since the LH aperture (ll.l"x22.3") is larger than the SH aperture (4.7"xll.3") and 
it is likely that some sources are so extended that unlike the LH slit, the SH slit could not 
cover the entire emission area, one sometimes need to multiply a scaling factor in order for 
the SH flux to agree with the flux of the LH part. The consequence of the aperture difference 
shows up clearly in IRAS 23239+5754, for which the SH spectrum must be multiplied by a 
factor of ~ 5.2 to agree with the LH flux at the overlap part and for the whole spectrum to be 
consistent with the previous MSX (Midcourse Space Experiment) and IRAS photometry. As 
shown in Figured! for IRAS 23239+5754 the SH slit covers less than half of the emission range 
while the LH slit covers most of the emission area, confirming a factor of ~ 5.2 difference. 
The flux matching between the SH part and the LH part is much better for IRAS 16456-3542, 
and 18354-0638. We just need a scale factor of ~ 1.2—1.3. 



See http://adsabs.harvard.edu/abs/2006sptz.prop30403J 
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The reduced spectra of these three sources are displayed in Figure El where the objects 
are ordered in terms of their evolutionary classes: post-AGB stars and PNe. Based on 
their ionized nebular emission lines, we classify IRAS 16456-3542 and 23239+5754 as PNe, 
consistent with the previous classifications in the literature. For the sources lack of nebular 
lines, they are commonly considered either as post-AGB stars or as extremely evolved AGB 
starsE In this context we classify IRAS 18354-0638 as a post-AGB star (see §3.1). 

The IRS wavelength range between 9.9 /im and 37.2 /im is full of spectral features arising 
from crystalline silicates. In order to clearly illustrate the crystalline silicate features, We try 
to subtract the underlying continuum. Instead of adopting a single-temperature blackbody 
spectrum, we simply fit the global continu um by a mathematical spline function. This 



approach was taken by lMolster et al.l (l2002bl ) in a comprehensive analysis of the ISO spectra 
of 17 oxygen-rich circumstellar dust shells surrounding evolved stars. For an overview of the 
IR spectral energy distribution (SED), previous photometric data, in the 2MASS J band, 
MSX A banc0 and IRAS 60 /im band, are supplemented. The continuum is then maximized 
and subtracted. In principle, the strength of the features might be underestimated since 
there is a possibility that very broad features are mistaken as continuum. After this global 
continuum subtraction, one more step is taken to remove the "local continuum" by a third- 
order polynomial. The so-called "local continuum" refers to the part that is above the global 
continuum while below the minimum of the features. This step would possibly further 
underestimate the strengths of the features. Figure [3] is an example for subtracting the 
continuum underneath the spectral features. 

Figure H] displays the continuum-subtracted spectral features of these three sources 
(IRAS 16456-3542, 18354-0638, and 23239+5754). Tabled lists the peak wavelengths and 
FWHM (full width half maximum) of the crystalline silicate features derived from a Gaus- 
sian fitting. The spec tral features are compared with the ISO spectrum of NGC 6302, a PN 



(IMolster et al.ll2002bl ). The features previously identified in the Molster sample are labeled 



by thin vertical lines. 

Probably due to the high resolution and sensitivity of Spitzer /IRS, several "new fea- 
tures" (labeled by thick vertical lines in Figure H]) appear to be present in the IRS spectra 
of these sources: (1) the 11.6 /im feature seen in IRAS 16456-3542, (2) the 14.4 /im feature 
(or the 14.6 /im feature) seen in IRAS 16456-3542 and IRAS 23239+5754, (3) the 15.4 /an 



6 The AGB and post-AGB distinction is made based on previous studies available in the literature as the 
Spitzer spectra provide no further decisive information for this classification. 

7 The MSX A (msxA), B x (msxBl), B 2 (msxB2), C (msxC), D (msxD ), and E (msxE) bands respectively 
peak at 8.28, 4.29, 4.35, 12.13, 14.65, and 21.34 /im jEgan fc Pricelll996l) . 
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feature seen in IRAS 23239+5754 and in IRAS 16456-3542^ and (4) the 25.47 /xm feature 
in IRAS 23239+5754 (or the 25.85 /xm feature in IRAS 16456-3542). These "new features" 
have not previously been reported in the literature. A more detailed study of a large sample 
of both astronomical and laboratory crystalline silicate spectra is necessary for confirming 
these "new" features. 



3. Results: Remarkably Rich Spectra of Crystalline Silicates at 10—37 /xm 

3.1. Overview of the Three Objects Rich in Crystalline Silicates 

The three sources which show abundant features of crystalline silicates are IRAS 16456- 
3542, 18354-0638, and 23239+5754. We discuss these objects below, in the order of their 
evolutionary stages. 



1. IRAS 18354-0638. The Spitzer /IRS spectrum of this source exhibits one nebular line 
(i.e. [SIII] at 33.5 /mi) which may suggest the start of a nebular phase. We hesitate to 
classify it as a PN because PNe usually exhibit > 3 nebular lines (while IRAS 18354- 
0638 displays only one emission line). We classify it post-AGB star@ In comparison 
with the other two objects (IRAS 23239+5754 and 16456-3542; see below) which display 
abundant crystalline silicate features and are classified as PNe, this object seems to be 
at an earlier stage of evolution (i.e. extreme AGB or post-AGB). 

2. IRAS 23239+5754- This object is a very young bipolar PN with a complex inner 
structure inside the bi polar lobes. The effective tempe rature of the central star is 
about T e fj w 40, 000 K (IPreite-Martinez fc Pottaschlll983l ). The mean temperature of 
the bulk dust in the circumstellar shell is ~ 200 K, as estimated from the overall IR 



spectral energy distribution (IZhang fc Kwoklll990l ). Several nebular lines (e.g., [SIV] at 
10.5 /xm, [Nell] at 12.8 /xm, [SIII] at 18.7 /xm, [SIII] at 33.48 /xm, and [Nelll] at 36.0 /xm) 
clearly show up in its IRS spectrum. This sour ce is well known for its pure fluorescent 
H 2 emission in the near IR wavelength region ( iRamsay et al.l 119931 ). It has long been 
suspected to be a binary. 



Also known as Hubble 12 at a distance of ~3kpc (jKingsburgh Barlowlll992l ). this 
PN, in particular its multiple coaxial rings in the whole-scale bipolar structure, has 



8 We canno t exclud e the possibility that this feature may be the same as the 15.2 /xm feature reported in 
Molster et all |2002bl ). 



3 We note that this star is not in the Torun catalog of post-AGB stars (jSzczerba et al.l 120071 ). 
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been well studied in the literature (see e.g. iKwok fc Hsial 120071 and references therein). 
The shape of this object reve als that it is probably at a very early PN phase, with a 
kinematic age of ~300 years (IMiranda fc Soh1ll989l ). The multiple rings extend about 
4"x 8". This size is comparable to the aperture of the IRS SH slit (3.7" x 11.3") and 
much smaller than the IRS LH slit (11.1" x 22.3"). As discussed in §2, this explains 
the discrepancy between the SH and LH fluxes. 

3. IRAS 16456-3542. This object is also known a s PN G347.4+05.8 and Hl-2. Its central 
star is a "Weak Emis sion-Line Star" (W ELS; iGesicki et al.l 120061 ) . With a dynamical 



age of ~ 1400 years (IGesicki et al. 120061 ) . it is a young PN (but more evolved than 
IRAS 23239+5754). IPhillipsI (120021 ) derived a dis tance o f ~2.4 7kpc, but la rger dis- 
tances have also been suggested (e.g., ~6.96kpc [Zhangl (119951 )]. ~9.01kpc [Phillips 
(120041 )] ). In addition to the crystalline silicate bands, this object also emits at the 11.3, 



12.7 and 16.4 /im PAH bands (see Figure [2]). 



3.2. Percentages of Crystalline Silicates 



To derive the mass percentage of crystalline silicates, we follow iMolster et al.l (l2002d ) 
to apply a simple dust emission model to fit the observed IRS spectra (including the con- 
tinuum underlying the crystalline silicate features). This model assumes that the dust shell 
is optically thin in the IR. The SED is then fitted with mass absorption coefficients K a bs(^) 
multiplied by blackbody functions B(T, v): 

F(i/)=£[fl(T 4 ,i/)x/& 1 ,(i/)xm,] . (1) 

The dust species considered here include: (1) amorphous silicate0 (2) forsterite, and (3) 
enstatite. For each species, we consider both a cool component and a warm component 
respectively at temperature T c and T w The mass of dust species i, m i} is a multiplication 
factor related to the total mass and only represents the relative mass ratio 1^1 



10 For IRAS 18354-0638 amorphous carbon fits better than amorphous silicates. For this object, the crys- 
talline fraction is the mass fraction of crystalline silicate relative to the total dust mass. 

11 The dust temperature is treated as a free parameter in the range of 40 K < T < 1000 K. We start with 
a single-temperature component and then add a component at another temperature to see if the reduced \ 2 
significantly drops. 

12 To derive the true dust mass, the knowledge of the distance of the object is required. The relative mass 
ratio of different components does not depend on the distance and can be correctly determined in this way. 
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The absorption coeffi cient of amorphous silicate is calculated from the optical constants 
of lOssenkopf et al.l (ll992f ) with Mie theory, assuming the dust to be spherical with radii 
a = 0.1 /im, the typical size of interstellar du S t@ For the minerals (forsterite and enstatite), 
the optical properties measured in laboratory depend on the size, shape , struc ture and 
chemical composition of the dust analog. We adopt those of iKoike et al.l (119991 ) because 
their coefficients closely resemble the IRS spectra of our sources. 

The relative mass of each component (amorphous silicate, forsterite and enstatite) is 
treated as a free paramet er. For enstat i te, we consider two sub-types: clino-enstatite and 
ortho-enstatite. Following iMolster et al.l (12002d ) . we assume equal mass of ortho- and clino- 
enstatite0 

We take the Levenberg-Marquart minimization algorithm to obtain the best-fit. The 
results of the SED fitting are shown in Figure and the parameters are listed in TableSJ The 
mass percentage of crystalline silicates is ~ 10 ± 2%, 82 ± 5% and 97 ± 2% respectively for 
IRAS 18354-0638 (see Footnote HDD , 23239+5754 and 16456-3542. Except for IRAS 18354- 
0638, the fraction of crysta lline silicate is much higher than the average of ~ 10% in evolved 
stars (jKemper et al.ll200ll ). and the amount of crystalline silicate even exceeds amorphous 
silicate. In IRAS 16456-3542, the fraction of crystalline silicate reaches ~ 97%, which exceeds 
the highest crystalline fractio n measured so far fr om the crystalline silicate emission spectra 
(~75% in IRAS 09425-6040 IMolster et al.l l20~0~ll ) f^l Even when we take into account the 
uncertainty (< 5%) in the modeling, the fraction of crystalline silicate in this source is one 
of the highest. We conclude that there must be some evolved stars in which the silicates are 
even completely crystallized. 

Under the assumption that the crystalline silicate in the evolved stars has only two forms: 
enstatite and forsterite, the mass ratio of the two forms m e /mf is approximately 4, 54 and 
24 respectively in these three sources, implying that the main form of crystalline silicate 
is enstatite i.e. MgSiC>3. Enstatite accounts for more than two-thirds of the crystalline 



13 The 9.7^m amorphous silicate band in our sources is better fitted in terms of small-sized grains than 
large ones. 



14 Enstatite emits at 20.65 /jm and 21.7/zm (jMolster et al.lll999l ). The relative strengths of these two 
features show a notable difference between clino-enstatite and ortho-enstatite: for clino-enstatite they are 
about equally narrow, while for ortho-enstatite the 21.7 (im feature is stronger and broader than the 20.65 /im 
feature. The 20.65 /im and 21.7 fim features are comparably narrow in the thr ee sources con s idered here and 
in AFGL 4106, the source in which these two features were first identified bv IMolster et al.l (|l999h . 

15 Several crystalline silicate absorption features were detected recently in a distant absorber at zw0.89 
toward PKS 1830-211, a gravitationally lensed quasar (|Aller et al.l 120121) . The crystalline fraction derived 
from the absorption spectrum of this object may be >95%. 
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silicate dust. This result is consistent with iMolster et al.l ( I2002cf ) who found the enstatite to 
forsterite mass ratio is between ~ 1 and ~ 11.4 for 12 evolved objects. 



Discussion 



The model fit presented in §3.21 shows that although most of the crystalline silicate 
features can be fitted reasonably well in terms of two silicate mineral types (i.e. forsterite 
and enstatite), some features are not well fitted (see Figure [5]). Most noticeably, the feature 
at 32.7/zm, prominent in all three sources and particularly strong in IRAS 18354-0638 (in 
which the 32.7/im feature is nearly as strong as the 33.5 /zm feature), is not accounted for by 
the current model consisting of forsterite and enstatite minerals. Similarly, the model fails in 
predicting the features at 29.6, 30.6 and 31.1 /u,m seen in all three sources (see Figure [5]). We 
note that the 29.6, 30.6 and 31.1 fim features seem to comprise a group as they have similar 
relative intensities in all three sources, with the 29.6 /im feature being the strongest and the 
31 .1 /xm feature b eing th e weakest. Such a phenomenon was also seen in the larger sample 
of IMolster et al.l ( j2002a[ ). One may speculate that these features arise from a dust mineral 
other than enstatite and forsterite. It is also interesting to note that the overall spectral 
features at 10-35 /im of IRAS 23239+5754 and IRAS 16456-3542 are similar to each other 
(probably because they are both very young PNe), except that the 29.6 fim feature is much 
stronger in the older PN IRAS 16456-3542 than that in the younger PN IRAS 23239+5754. 

In summary, the modeling approach discussed in §3.21 assumes three grain types: (1) 
spherical amorp h ous s ilicates of radii a = 0.1/im with the indices of refraction taken from 
Ossenkopf et al.l (119921 ); (2) crystalline forst erite and (3) cryst alline enstatite for which the 
mass absorption coefficients are taken from Koike et al. ( 19991 ) 1^1 For each grain type, two 
temperatures are assumed (i.e. a cold component and a warm component) unless a single 
temperature model results in more or less the same x 2 . To examine whether the silicate 
crystallinity derived in §3.21 is dependent on the modeling t echnique, w e take a diff erent 
approach: we follow the approach documented in detail by iLisse et al.l (120071 120081 ) and 
incorporate more sizes and more silicate mineral compositions. The detailed mineralogy of 
the dust is determined by comparing the observed Spitzer/IRS spectrum to a wide range 
of plausible silicate minerals with high-quality mid-IR transmission and thermal emission 
laboratory spectra. The calculated model flux from each mineralogical component is the 
product of a blackbody of characteristic temperature times the experimentally determined 
absorption efficiency for each sized particle, summed over all particles in the dust population. 



16 For crystalline forsterite and enstatite we do not need to assume grain size, shape, etc. 
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For illustration, in Figure |6] we show the Spitzer/IRS spectrum of IRAS 16456-3542 and the 
model spectrum as well as the appropriately scaled contributions from each grain species. 
The overall fit is reasonably good, although not as good as that described in §3.21 (see 
Figure E]). The derived silicate crystallinity is ~ 90%, consistent with the silicate crystallinity 
of ~ 97% derived in MM 



Finally, we note that various techn iques have been proposed in the literature to fit the 
crystalline silicate emission spectra (e.g.jMolster et al. 1 12002 a; Ho nda et a ISargent et al 
20061 ; iLisse et al.ll2007t iGielen et al.ll2008l ). All these approaches treat d ust size, comp osition 
and temperature as free parameters. This is physically incorrect (see iLi et al.l 120041 ) . The 
equilibrium temperature of dust is sensitive to its size and composition: smaller grains have 
higher temperatures, and silicates with a higher Fe/Mg ratio have higher temperatures. The 
dust temperature is also sensitive to the distance from the dust to the illuminating source 



sec 



Li et al.l 120041 ) . Therefore, one cannot treat the dust temperature as an adjustable pa- 



rameter independent of the dust composition and size. The dust mass is very sensitive to 
the dust temperature. The crystalline fraction derived from all the approaches in which the 
dust temperature is treated as a free parameter should be taken with caution. The correct 
approach should physically calculate the dust temperature for a give n dust size of a given 
composition at a given distance from the illuminating source (e.g., see 
Chen et al.ll2007h . 



Li fc Greenberdll998 



5. Summary 

We report the 9.9-37.2 /zm spectra of three evolved objects (post-AGB stars and PNe) 
using the high-resolution mode of the IRS instrument on board Spitzer. These objects exhibit 
rich sets of crystalline silicate emission features. These features are fitted by a mixture of 
three dust species: amorphous silicate, crystalline forsterite, and crystalline enstatite, with 
each species consisting of a warm component and a cold component. It is found that the 
crystalline silicates are dominated by crystalline enstatite. The crystalline silicate spectral 
modeling shows that IRAS 16456-3542, has a degree of silicate crystallinity (i.e. mass fraction 
of silicate dust in crystalline form) of ~ 97%, the highest to date in crystalline silicate sources. 
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Table 1. Infrared properties 



IRAS Name 


Type 


2MASS (magnitude) 




MSX (Jy) 






IRAS (Jy) 








J 


H 


K 


A 


c 


D 


E 


12 (im 


25 )im 


60 /an 


18354-0638 
23239+5754 
16456-3542 


Post-AGB 
Young PN 
PN 


8.55 
10.23 
11.83 


5.85 
9.83 
11.48 


4.51 
8.81 
10.53 


9.10 
10.99 
2.26 


19.60 
16.95 
3.84 


24.40 
18.53 
7.06 


49.10 
50.81 
22.19 


18.90 
20.86 
4.88 


68.00 
71.22 
29.61 


39.20 
35.79 
10.96 
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Table 2. Journal of the Spitzer /IRS observations 



Star No. IRAS Name MSX6C Name Date observed Data Version Comment 



1 18354-0638 G025. 5470-00.0529 04/20/07 S16.1.0 

2 23239+5754 Glll.8777-02.8509 09/09/06 partial S16.1.0 poor LH nod2 

3 16456-3542 G355. 7334-03.4721 03/25/07 partial S16.1.0 




Fig. 1. — Spitzer /IRS SH (small rectangle) and LH (large rectangle) slits superimposed on the 2.159 /im 
2MASS K band image (bright black and white) and the 21.34 /im MSX E band image (faint extended). 
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Table 3. Inventory of crystalline silicate features of IRAS 16456-3542, 18354-0638, and 

23239+5754 



IRAS 18354-0638 | IRAS 23239+5754 | IRAS 16456-3542 
A peak FWHM A pcak FWHM A pcak FWHM 



19 


45 





83 


20 


65 





34 


21 


66 





47 


23 


21 


1 


21 


24 


49 





23 


27 


73 





85 


29 


19 





85 


30 


54 





29 


31 


11 





19 


32 


78 





95 


33 


57 





88 


35 


41 


1 


36 



10.63 0.11 

11.26 0.39 

13.54 0.22 

13.75 0.25 

14.18 0.12 

14.40 0.13 

14.63 0.08 
15.45 0.33 
16.10 0.60 
17.47 0.99 
18.00 0.53 

19.59 0.54 
20.70 0.31 

21.60 0.45 
23.09 0.54 
23.74 0.49 

24.58 0.39 
25.47 0.55 
27.72 1.30 
29.53 0.86 

30.59 0.36 
31.12 0.14 
32.81 0.50 
33.80 0.93 

35.64 1.55 



10.60 0.24 
11.22 0.17 

11.61 0.27 

13.55 0.25 

13.75 0.23 

14.39 0.26 

14.61 0.19 

15.39 0.31 

16.10 0.30 

17.43 0.39 

18.96 0.42 

19.61 0.50 

20.65 0.29 

21.60 0.53 

23.07 0.50 

23.68 0.53 

24.52 0.24 

25.85 0.30 

27.88 1.14 

29.58 0.66 

30.57 0.27 

31.13 0.14 

32.80 0.42 

33.72 0.62 

35.55 1.14 



Note. — Compared w ith the crystalline silicate feature list 
of iMolster et al.l (|2002al ). t he features in Roman are at the 
same wavelength as that of IMolster et al.l (|2002af) . The fea- 
tures peaking at appreciably different wavelengths are listed 
in italic. Those listed in bold face are new features, not iden- 
tified previously. 
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Fig. 2.— Spitzer/WS spectra of IRAS 18354-0638, 23239+5754and 16456-3542 (in the order 
of evolutionary stage from post-AGB phase to young PN phase) which are rich in crystalline 
silicates. 
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Table 4. Dust parameters derived from crystalline silicate spectral modeling (see §3.2[) 



Source 


T cf 




T 

± ce 


T T 

- 1 we - 1 ca 


T 


m e /rri{ 


Crystallinity(%) 


18354-0638 


79 




87 


••• 85 


213 


4 


10 ±2 


23239+5754 




129 


81 




150 


54 


82 ±5 


16456-3542 


67 


125 


72 


144 


158 


24 


97 ±2 



Note. — From left to right: the source, the temperatures of the cool forsterite component 
(Tgf), the warm forsterite component (T wf ), the cool enstatite component (T ce ), the warm 
enstatite component (T we ), the cool amorphous silicate component (T ca ), and the warm 
amorphous silicate component (T wa ), the mass ratio of enstatite to forsterite (m e /mf), and 
the derived silicate crystallinity (i.e. the mass fraction of silicate dust in crystalline form). 




10 20 30 40 50 60 

Wavelength (| im) 

Fig. 3. — Example of continuum subtraction. The global continuum is fitted by eye-guided spline simulation 
(dash line) of the observed spectrum (solid line), with the addition of three photometry data in the 2MASS 
J band, MSX A band and IRAS 60 /im band denoted by red triangles. The local continuum is fitted by a 
3-order polynomial with the continuum-subtracted spectral features shown in the inset. 
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Fig. 4. — Rich crystalline silicate fea tures in IRAS 16456-3542, 18354-0638, and 23239+5754. The features 
already listed in lMolster et al. ( 2002a ) are labeled with thin solid black line. The "new features" not reported 
in the literature are labeled with t hick so lid red line. Those with noticeable shifts of peak wavelength 
(compared to that of iMolster et al.l (|2002al) ) are labeled with dotted red line. 
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Fig. 5. — Comparison of the model spectra (red solid line) with the Spitzer/IRS spectra (gray 
solid line). Three types of silicates are considered: forsterite (blue dashed line), enstatite 
(ortho-enstatite and clino-enstatite; green dot-dashed line), and amorphous silicate (red 
dotted line). For IRAS 18354-0638 we consider amorphous carbon instead of amorphous 
silicate. 
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Fig. 6. — Modeling the Spitzer/IRS spectrum of IRAS 16456-3542 with the approach of 



Lisse et al.l ( 120071 . [20081) which incorporates a power-law size distribution of ~ a with 



0.1 (Um< a < 10 /im for the grain size and various crystalline and amorphous olivine and 
pyroxene grain materials. The solid lines beneath the IRS spectrum show three components: 
the upper red and green lines for amorphous carbon and carbonates, the middle cyan lines 
for pyroxenes and sulfides, the lower blue lines for olivines. 



